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Abstract 

We tune the calculation of central exclusive Higgs production to the recent CDF central exclusive 
dijet data, and predict the cross section for the exclusive production of Higgs boson at the LHC. It 
is always below 1 fb, and below 0.3 fb after experimental cuts. 
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1 Introduction 



The possibility of producing a Higgs boson in a diffractive event has been suggested a long time ago 
Pfl E] • The exclusive process was then recognized as a gold-plated signal for a light Higgs boson [3j , 
as it would show up irrespectively of its decay products. Many theorists built models predicting 
the rate of production, and because there is no factorisation theorem for production in a gap, the 
estimates varied wildly. Recently, CDF measured diffractive production of jets in a mass range up 
to 130 GeV, so that one can now identify the ingredients needed to describe the data. 

The first attempt to embed Higgs boson production into a pomeron appeared in [2], which used 
specific non-perturbative ingredients to model the pomeron. The first perturbative calculation was 
performed in [3], at the price of the introduction of an unknown proton impact factor. Potentially 
large screening corrections were evaluated in [5] in an eikonal formalism, and have been extended to 
any unitarisation scheme in (6J [7] . Finally, the necessity of a Sudakov form factor was noted in [8] . 

All these ingredients were incorporated in a series of papers by the KMR group [9], who used 
specific prescriptions for the Sudakov form factor, for the screening corrections, and for the impact 
factor. In a previous paper [10] devoted to dijet exclusive production, we have shown that their 
choices are not the only ones possible, and we have proposed alternatives to several of them: 




• Partonic amplitude: although the calculation lies partly in the infrared region, one usually 
calculates the amplitude in perturbative QCD by computing the exchange of two gluons between 
two quarks, Fig. [l](a). The second (screening) gluon is present in order to prevent a colour 
flow in the t channel and allows the production of a colour singlet in the final state. In the 
present calculation the momentum lost by the proton is not neglected with respect to the loop 
momentum k as in [TO], and we have shown that this approximation overestimates the cross 
section by a factor 2. An exact transverse kinematics is thus kept. 

• Impact factor: the infrared divergences of the QCD amplitude disappear because the protons 
is colour- neutral. This introduces a proton impact factor, as illustrated in Fig. [TJ(b). In the 
present paper,it is parametrised by an off-diagonal skewed unintegrated gluon density [TOj, [JT] 
that allows for the exchange of soft gluons in the t channel. 

• Sudakov form factor: the virtual corrections, coming from the transition of gluons from low 
transverse momenta in the proton to hard ones in the Higgs-boson vertex, are accounted for 
by this factor, as shown in Fig. [TJc. The vertex corrections give double and single logarithms 
which must be resummed and the main uncertainty in the calculation comes from the constant 
terms that cannot be resummed although they can be large. 

• Gap survival probability: additional pomeron exchanges between the initial and final protons 
can occur and this leads to screening corrections [5], as shown in Fig. fTl (d) 



1 



We also showed that the calculation still largely lies in the non-perturbative region, so that more exotic 
choices are also possible |12| . The uncertainties coming from these various ingredients were studied, 
and in the following we shall use the CDF dijet results [10] to reduce the theoretical uncertainties on 
the Higgs-boson exclusive cross section. 

The first section of this letter summarizes the calculation of the Higgs-boson exclusive cross 
section in the light of the dijet exclusive cross section. The second section presents the results using 
our new knowledge of the different ingredients. 



2 The Higgs-boson exclusive cross section 
2.1 Kinematics and partonic amplitude 

The basis of the calculation is the partonic subprocess qq — >■ q + H + q, and involves the evaluation of 
the diagram shown in Fig. [2] and the corresponding cross diagram as well as the diagrams in which 
the Higgs is emitted from the other gluon. We neglect the quark masses, their initial transverse 




Figure 2: Kinematical conventions for the Higgs exclusive production at the partons level. 



momenta and write their four-momenta as q^ and p^ with s = 2p ■ q. The total momentum transfers 
to the initial quarks are —k\ and k% respectively and the momenta are decomposed onto the incoming 
quark directions using a Sudakov decomposition: 

& = ap^ + (3q^ + k* 4 , 

k*t = a lP » + p x <f + k£, (1) 
K = a 3 p» + foq» + k£, 

where all transverse vectors are denoted by bold letters. The produced Higgs boson momentum is 

= (ai - a 3 )p M + (ft - /W + (ki - ks)". ( 2 ) 

The kinematical region of interest is multi-Regge kinematics where longitudinal momentum losses 
are small and obey 

|k 2 | |k 2 | 
1 » ai,pi > Ljkl , * = 1,3; a,/3~J — [ . (3) 
s s 

Let us also define the momenta of the colliding gluons 

= (* + h)", 4 = {k + k 3 y. (4) 
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At lowest order, the amplitude is purely imaginary and can be calculated using cutting rules. 
The parton-level amplitude then becomea^J in the high-s limit [lj: 

where g s is the coupling constant between quarks and gluons, which will be later re-absorbed in the 
impact factor, C=2/9 is the colour factor for the two-gluons exchange and V^ v is the Higgs-gluon- 
gluon effective vertex. The effective coupling between gluons and the Higgs boson can be obtained 
by calculating the diagrams in Fig. [31 This was done in [3j where it was shown that the vertex can 




Figure 3: The Higgs-gluon-gluon vertex. Because the coupling of the Higgs to quarks is proportional 
to the quark mass, the top quark contribution dominates. 



be written 
where 
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with g the hard coupling constant, Gf the Fermi constant and m/j the Higgs mass. The expression of 
N\ and A^2 come from the integration over the quark loop in the particular kinematics of production 
in a rapidity gap k|,k§ <C m%- Defining a = m 2 H /m 2 , one obtains 



1+1 



arctan 



The second part of the vertex, that depends of N2, can be neglected because the tensor structure will 
always give contributions proportional to k\ or k 2 that are small compared to the factors s generated 
by the first tensor structure, and because |AT 2 | 2 is always less than 30% of |iVi| 2 for Higgs masses 
below 1 TeV. 

Considering now the phase space available for the final Higgs boson and the above expression for 
the vertex, one can write the partonic differential cross section in the large-s limit as 



2^8 
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(9) 



The final delta function is used to put the Higgs on shell and to perform the integral over 0:3. 
One can show using Eq. Q and the on-shell condition on the final quarks that 



777 



a 3 



H 



s Pi 



(10) 



The above expression, Eq. ([9]), is singular at the partonic level when one of the exchanged gluons 
goes on-shell, i.e if \t 2 K 2 K 2 — > 0. This divergence disappears when one takes into account the fact 
that protons are colour singlets. This introduces an impact factor which furthermore describes the 
response of the proton to a transfer of momentum. 



4 Note that ref. [4] is missing a factor 2 in the expression of the amplitude. 
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2.2 Soft and Higher-Order Corrections 



The lowest-order partonic amplitude is the only piece of the calculation that can be exactly calculated. 
The higher-order corrections (Sudakov form factor) and the soft corrections (impact factors and 
screening) can only be estimated. One can also use other exclusive data to constrain the theory, 
and the only other input in a similar kinematical range is the dijet cross section from the TeVatron. 
Hence we shall use our previous study [TO] to constrain the calculation. 



2.2.1 Impact Factor 

The wave- function overlap between the initial-state and final-state protons, after an exchange of 
i-channel gluons of transverse momenta k a and k^ and of longitudinal momentum fractions x a and 
Xb, leads to an impact factor $(x a , Xf,; k a , k&) that has the properties [13] 



$(x , x b ; 0, kj,) = $(x a , x b ; k a , 0) = 0. (11) 
The amplitude ^ for Higgs-boson production becomes for proton-proton scattering 

• Mpp = ff4 2(2^ / / ^p( 4 9%)^ l/ ( 4 ^) $ (/3 ) /3i;k,Ki)$(a,a 3 ;k,K 3 ) (12) 

To this, one must add an s-dependent factor that comes from the gluon content of the proton 
(or equivalently from the pomeron trajectory), so that $ is a function of the longitudinal momentum 
fractions of the gluons. As the spectator gluon of momentum k has a smaller momentum fraction 
than that of the gluon linked to the Higgs vertex, we assume that the latter momentum fraction is 



dominant, and take a = f3 = in Eq. (12). 

Although the property (11) is true in general, the precise form of $ is unknown. In [10], we 
tested several possibilities, and found that the shape of the Et distribution of the CDF dijet data 
was best reproduced if we assumed an x-dependence similar to that coming from unintegrated gluon 
distributions. Hence we built the impact factor using three ingredients: 

• The diagonal unintegrated gluon distributions J-(x,Q 2 ) is built in [llj. They are built from 
a hard part, based on differentiation of conventional gluon densities (GRV [13], MRS [TO] or 
CTEQ [16] LO fits), and from a soft part, based on soft colour-singlet exchange inspired by the 
dipole form factor and allowing partons to enter the IR region. Four variants of the soft part 
were used, and their parameters were fixed by fitting the proton structure function i*2 P to the 
HERA data in a region of virtuality < Q 2 < 35 GeV 2 and for x smaller than 10 -2 . They are 
respectively named FIT-1, FIT-2, FIT-3 and FIT-4 [171 ITS], 

• The skewed unintegrated distributions are obtained from the above by assuming that skewness, 
i.e. the non-zero longitudinal momentum transfer, can be taken into account by using an 
effective longitudinal momentum fraction x g = 0.41 Xj |17|. 119]. 

• The fully off-diagonal impact factors are obtained by multiplying the above by a factor repro- 



ducing Eq. (11) and going to 1 in the diagonal case, and by an exponential factor reproducing 
the proton t slope and shrinkage. One also chooses the scale Q to be the average of the incoming 
transverse momenta Q 2 = ( ka + kj > ) 2 . 



One then obtains for the upper proton: 



HO, ft; k, = -F (0.41/3,, (k + i kl ) 2 ) 6 " J 



(13) 
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with B =A GeV" 2 , a'=0.25 GeV" 2 and x = 3.4 10" 4 [TO]. 

These parametrisations provide a very good description of the charm contribution to the proton 
structure function i 7 ^, of the longitudinal structure function and of diffractive vector meson 
production [171 [20] as well as of the dijet exclusive cross section at the TeVatron [TO] . 



2.2.2 The Sudakov Form Factor 



During the exchange of gluons producing the Higgs boson, the gauge field goes from a long-distance 
configuration to a short-distance one. This is the situation in which, for running as, large logarithms 
arise from virtual diagrams such as those of Fig. [TJ (c) . These logarithms become infinite if the external 
legs go on-shell, and it is known that these infinities are canceled by the bremstrahlung diagrams, 
which are easier to calculate. Hence, in order to evaluate the vertex corrections to g*g* —> H, one 
can use the standard formula |21j . assuming a point-like ggH vertex: 



2 n2\ 



-S{f?,l 2 ) 



, 5(/i 



2 o2\ 
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where A = |q|//x and P gg and P gq are the splitting functions 
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The amplitude becomes 
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(14) 
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The Sudakov form factor is a loop correction to the vertex and includes an integral over the mo- 
mentum in the loop, i.e. over both the fraction of light-cone momentum z and the transverse 
momentum q 2 . The integration builds up double logarithms in the kinematical region where the 
virtual emission is soft and runs from a soft to a hard scale, respectively, £ 2 and /i 2 . The upper scale 
is chosen according to the standard discussion of the Sudakov form factor where the hard point-like 
vertex is characterized by only one scale, the mass of the centrally produced system mx ■ This is the 
situation in the Higgs boson production for which the log structure of the Sudakov form factor can be 
calculated. It has been estimated to single-log accuracy in [22] and recently in [23], the upper limit 
of the z-integration is then of the order of ran- The lower scale £ 2 is proportional to the virtuality 
from which the evolution starts ftj and this leaves an uncertainty related to the constant terms which 
cannot be obtained from bremstrahlung, and which are known not to exponentiate. For the purpose 
of this study, we perform the z integration in the argument of the Sudakov correction S(/j?,£ 2 ) and 
write its structure in three terms. After the second integration considering a running as, the first 
will give the log contribution, the second the (log) (log) contribution and the third a constant term: 



with 
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Physical Quantity 


Allowed Values 


Fit to the 


1 a Values 






CDF data 




Parametrisation of the 


[321123 


FIT-4 


FIT-4 


unintegrated gluon distribution 








Constant terms in Sudakov exponentiation 


C/2 to 2 C 


C=1C 


C/2 to 2C 


Gap survival probability 


4.5-18% 


10 % 


10% 



Table 1: Allowed range of parameters, and parameters reproducing the central values of the CDF 
exclusive dijet data. These will be used to predict the Higgs-boson exclusive cross section. 



We have kept the latter contribution as in other exclusive models but we allow its variation by a 
factor 2 up or down. This constitutes the uncertainty coming from the constant terms in the Sudakov 
form factor. The effect of this change is an overall factor. 

Note that perturbation theory suggests that the Sudakov form factor should be included in the 
definition of the unintegrated gluon distribution. We do not do this here, because the use of an 
unintegrated gluon distribution is only an educated guess leading to the impact factor. There is 
no reason to believe that the two objects are identical, as the exchanged spectator gluon cannot be 
assigned to only one of the initial protons. 

2.2.3 Gap-Survival Probability 

There is no factorisation theorem in diffractive exclusive production, so it may be important to take 
initial- and final-state exchanges into account. This is at best tentative, as one really needs to model 
the multiple-pomeron exchanges which account for the purely elastic cross section, i.e. one needs to 
overcome the difficulty that lead the analytic S-matrix theory to a standstill. 

Most of the present estimates are based on an eikonal scheme of unitarisation, following the 
arguments presented by Bjorken a while ago They have been recently generalised to an arbitrary 
unitarisation scheme [7J, for the case in which the production process happens at short distance. 
The gap survival probability can be calculated as the S"-matrix element, which factors in impact- 
parameter space and depends on the variables conjugate to k\, k^ and k\ + k^. Such construct, and 
its transform back to momentum space, has never been implemented, and one usually treats the gap- 
survival probability as an overall factor. Several estimates can be found in the literature [24J, from 
5% to 15% at the TeVatron energy, and about a factor 2 lower at the LHC. These are conservative 
estimates, as the gap survival at the LHC could be as large as 25 % or smaller than 1 % [7]. In the 
present study and in Sec. |4[ we assume a gap survival probability at the LHC two times smaller that 
the one used at the TeVatron as suggested by eikonal models [24]. 

2.3 First Estimation of the Higgs Boson Exclusive Cross Section 

Putting together all the above ingredients, we can estimate the Higgs boson cross section for sets 
of reasonable theoretical parameters given in the second column of Table [TJ and we see that 
modest changes in the theoretical input can produce large differences in the predicted cross section. 
The Higgs boson exclusive cross section in this model is below one femtobarns at the LHCj^J The 
uncertainty is of one order of magnitude and the slope as a function of the Higgs mass is steeper 
than in the existing similar models [J2 [25] . 

Similarly small cross section has been obtained in references |25l I26| . 
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Figure 4: Prediction of the Higgs exclusive production cross section at 14 TeV without any experi- 
mental cuts and without tuning to CDF dijet data. 



3 From Dijet to Higgs 

One can in fact do a somewhat better job by constraining the model to reproduce existing data. 
The process of interest is the dijet exclusive production for which one of the diagrams is shown in 
Fig. plb and for which the CDF collaboration published a measurement in 2008 [30]. These are 




(a) (b) 

Figure 5: (a) Higgs exclusive production and (b) Dijet exclusive production 

important data for two reasons, firstly the calculations are expected to be almost identical, i.e. they 
differ only because of the hard vertices that produce the final particles. Indeed, diagrams in which 
the final-state gluons are produced by two different t-channel gluons can be safely neglected in the 
present kinematics as it contains one more hard propagator. In [27] , we have evaluated its maximum 
contribution and have showed that it is safe to consider only those of Fig. [5] Secondly, the dijet 
system reaches masses of the order of 130 GeV, i.e. the region of mass where the Higgs boson is 
expected. The CDF data thus provide a good opportunity to tune the calculation of quasi-elastic 
production and to reduce the uncertainties that affect this family of process. The model of production 
of exclusive dijet is the same as the one presented in the first section and developed in |10j . However, 
the sources of uncertainties on both processes are different due to the structure of the Sudakov form 
factor. One knows that a typical double-log-enhanced correction gets sizable logs only for a point-like 
vertex, or in other words, only if there is a large virtuality inside the effective vertex. This is the 
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case in the Higgs boson production where the virtuality of the top quarks in the vertex shown in 
Fig. [3] is suffisuantly large in front of q 2 . On the contrary, in the dijet case, one can end up with 
the two situations pictured in Fig. [6] In the left picture, the flow of transverse momentum q 2 <C k 2 . 




Figure 6: The two different kinematical regimes for colour-singlet dijet production. Left: The effective 
point-like vertex q 2 < k^. Right: The resolved vertex q 2 ~ k|. 



in the loop does not change the amplitude while in the right one, the momentum q 2 ~ k 2 is enough 
to resolve the vertex. The transverse logarithm does not build up leading to a complete change of 
the log structure of the amplitude. Consequently in the dijet exclusive production, the upper scale 
must be of the order of k 2 , as it marks the limit beyond which the double-logs do not build up. This 
raises the question whether the dijet exclusive production data can be used to tune the Higgs boson 
calculation due to the extra uncertainty coming from the Sudakov form factor. 

One can note that the gap survival probability is an overall factor that multiplies the cross section 
and a source of uncertainty similar to the one coming from the constant terms of the Sudakov form 
factor. Both uncertainties can be considered as a single one and cannot be reduced independently 
and this disavantage is reduced if one considers the ratio between the Higgs boson and the dijet 
exclusive cross section. It removes the uncertainty coming from the dependence of the cross section 
on the energy behaviour of the gap survival probability and underlines the importance of dijet data 
at the LHC. The uncertainty on the constant terms is reduced but not removed as the hard scale 
used in both processes is different (fj, ~ Et in the dijet case and fi ~ run in the Higgs case). This 
means that, whatever is the precision reached by the dijet data at the LHC, the precision on the 
theoretical prediction for the Higgs boson exclusive cross section will be smaller. 

The comparison between different parametrisations of the dijet exclusive production model and 
the CDF data have led to a selection of parameters that reproduce the measurement well, as shown 
in Fig. [7} The central curve corresponds to one of the different choices of parameters that reproduce 
the central values of the data points and the band corresponds to parametrisations that are within 
the 1 a errors, also given in Table [TJ Using the knowledge of the TeVatron dijet data, we now 
keep the parametrisation of the unintegrated gluon distribution, the constant terms in the Sudakov 
exponentiation and the gap survival probability, we introduce them in the Higgs boson exclusive 
production calculation and present a prediction for the cross section using the parameters from the 
third and fourth columns of Table [TJ 

4 Final results 

We now present the results for the cross section at the Higgs boson level but experimentally, what 
is detected in the central detector is the decay products of the Higgs, a bb pair in the mass range 
considered. To obtain the cross section, the final jets have a priori to be reconstructed via a jet 
searching algorithm and this should be taken into account as the transverse energy of the parton 
is not equal to the transverse energy of the jet due to the energy lost when reconstructing the jet. 
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Figure 7: CDF Run II data on exclusive dijet, and a band of parametrisations that all go through the 
data. The curve corresponds to a choice of parameters that best fit the central values and the band 
corresponds to parametrisations that go through the 1 a errors. Both include the effect of splash-out 
discussed in [10J. 



The effects of the so-called splash-out were studied into details for gluon production in [18] but are 
not included here in the Higgs case. They are of course included in Fig. [7} Indeed, one can use the 
very forward detectors to reconstruct the mass of the central system from the measurement of the 
fraction of momentum lost by the initial hadrons. One in fact scans the mass domain to obtain the 
cross section without using the information coming from the central detectoij^J 

At the TeVatron (- v /s=1960 GeV) for the experimental cuts of the first column of Table [2j the 
Higgs boson exclusive cross section is tiny, of the order of 2 attobarn, for m#=120 GeV and the 
parametrisation of the central curve. That makes the Higgs boson impossible to observe at the 
TeVatron, at least if it is produced in an exclusive reaction. For the LHC, we present in Fig. |8]three 



Cuts 


TeVatron [30J 


LHC[31 


LHC [28 




[0.03,0.08] 


[0.002,0.2] 


[0.005,0.018] 


Pi 




[0.002,0.2] 


[0.004,0.014] 


\yx\ 


[-2.5,2.5] 


< 1 


< 0.06 


Gap size 


> 3.6 






M x 




> 50 GeV 


[80GeV,160GeV] 



Table 2: Experimental cuts for the TeVatron and the LHC. 



sets of curves at 7, 10 and 14 TeV. The first row corresponds to the cross section without experimental 

6 Actually, one has to match the presence of hadrons in the forward detectors with an Higgs boson event in the 
central detector but the mass of the centrally-produced system and the cross section can be obtained using forward 
detectors only [25]. The effect of the Higgs decay and the reconstruction of the bb final state were studied in the same 
model in reference [29] . 
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Figure 8: Higgs-boson exclusive cross section at 14 TeV (left), 10 TeV (center) and 7 TeV (right). 
The first row gives the results without cuts. The second row corresponds to the cuts of the third 
column of Table [2] and the third row to those of the fourth column 



cuts, the second for cuts resulting from the use of the forward detectors FP420+RP220 [31], given in 
Table 2 and the third for a more restrictive set of cuts [28J used in Monte-Carlo studies, also given in 
Table 2 The FP420 detector would be located in a small space left on the beam line of the LHC and 
be able to collect protons flying at 10 to 5 mm from the beam. With the addition of RP220 it should 
be possible to measure precisely the momentum of protons that have lost a fraction 0.002< £ <0.2 
of their initial momentum with an acceptance close to 100%. This measurement would determine 
the mass of the centrally produced system and allow a clear identification of the exclusive events by 
associating the information from both central and forward detectors. 

In the left column that gives the cross section without any experimental cuts, one can see that 
the uncertainty was reduced by a factor 5 compared to Fig. [4] by using the CDF data. The effect of 
the cuts is important and cannot be neglected: the application of the FP420+RP220 cuts reduces 
the cross section by a factor three. The rapid fall down of the cross section in the bottom row is 
also purely due to the effect of cuts. If one considers only the use of FP420, there is a cut-off on the 
longitudinal momenta lost by the proton, 

0.002 < |a 3 |,/?i < 0.02. 
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With 1 03) roughly equal to m 2 H /sfi\ if one considers a Higgs boson with a mass of 200 GeV in a 
process with y/s = 10 TeV, the only possible values for | CV3 1 stand between 0.02 and 0.2, and are 
outside the range of the forward detector. For collision of protons below 14 TeV, the use of the roman 
pots of RP220 in addition with the very forward FP420 is needed to increase the acceptance at large 
masses. This effect disappears at higher energies. 

5 Conclusion 

In a previous paper |10j . we have proposed a new model for quasi-elastic production that is similar 
to the KMR model, but differs from it in the implementation of the different parts of the calculation: 
we use an impact factor with the correct IR properties, exact transverse kinematics and include a 
parametrisation of the nonperturbative region. Our model was used here to obtain the cross section 
of exclusive Higgs-boson production at the LHC and the corresponding uncertainties. The vertex 
corrections, i.e the Sudakov form factor, was one of the main sources of uncertainty in the dijet case 
but in the Higgs-boson production case, a recent NLO calculation sets its structure and scales. The 
other pieces of the model, the parametrisation of the unintegrated gluon density, the gap survival 
probability and the constant term in the Sudakov form factor were chosen to fit the CDF dijet 
data. Using these results, the Higgs-boson quasi-elastic cross section is found to be between 0.07 and 
0.2 femtobarns for a Higgs boson of mass 120 GeV and \/s = 10 TeV considering the kinematical 
cuts of RP220 and FP420. This is a factor 2 lower than the KMR estimate [9]. 

Finally, we want to stress that the theory of exclusive production is not under full theoreti- 
cal control. Indeed, exclusive production does not factorize, and wanders appreciably in the non- 
perturbative region. Hence, many of the assumptions made here are rather conservative as they 
derive from perturbative QCD. The constraint that leads to a reasonable agreement between models 
is the reproduction of the CDF dijet data. 
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